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ABSTRACT The interpretation electrical phenomena in biomembranes is usually based on the assumption that the 
experimentally found discrete ion conduction events are due to a particular class of proteins called ion channels while the 
lipid membrane is considered being an inert electrical insulator. The particular protein structure is thought to be related 
to ion specificity, specific recognition of drugs by receptors and to macroscopic phenomena as nerve pulse propagation. 
However, lipid membranes in their chain melting regime are known to be highly permeable to ions, water and small 
molecules, and are therefore not always inert. In voltage-clamp experiments one finds quantized conduction events 
through protein-free membranes in their melting regime similar to or even undistinguishable from those attributed 
to proteins. This constitutes a conceptual problem for the interpretation of electrophysiological data obtained from 
biological membrane preparations. 

Here, we review the experimental evidence for lipid ion channels, their properties and the physical chemistry underlying 
their creation. We introduce into the thermodynamic theory of membrane fluctuations from which the lipid channels 
originate. Furthermore, we demonstrate how the appearance of lipid channels can be influenced by the alteration of 
the thermodynamic variables (temperature, pressure, tension, chemical potentials) in a coherent description that is free 
of parameters. This description leads to pores that display dwell times closely coupled to the fluctuation lifetime via 
the fluctuation-dissipation theorem. Drugs as anesthetics and neurotransmitters are shown to influence the channel 
likelihood and their lifetimes in a predictable manner. We also discuss the role of proteins in influencing the likelihood of 
lipid channel formation. 
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1 Introduction 

Biological membranes mainly consist of a lipid bilayer into 
which proteins are embedded. The mass (or volume) ratio 
between proteins and lipids ranges between 0.25 (lung sur- 
factant) and 4 (purple membrane of halobacteria). Typical, 
biomembranes display a protein-lipid ratio of approximately 
one. This includes the extra-membraneous parts of the proteins 
and surface-associated protein. Therefore, the intra-membrane 
parts of the proteins represent a smaller fraction of the central 
part of the membrane. Thus, even in densely crowded biologi- 
cal membranes most of the in-plane membrane area consists of 
lipids. 

Textbook models about the electrical transport properties of 
membranes consider the lipid bilayer as an electrical insula- 
tor. For instance, the Hodgkin- Huxley model of the nervous 
impulse considers the membrane as an inert capacitor and 
attributes currents to conductance changes occurring in ion 
selective protein channels 1 1 1. The activity of protein ion chan- 
nels is associated to quantized current events of order of 1-50 
pA (pico Amperes) at a voltage of order lOOmV corresponding 
to channel conductances of about 10-500 pS (pico Siemens). 
The typical open-dwell time of such channels is of the order of 
a few milliseconds. 

Interestingly, however, synthetic lipid membranes close to 
phase transitions are not inert but are very permeable for small 
molecules [2], ions ||3] |4] |5] 13 and water [6j. Further, they 
display stepwise conductance events very similar to those of 



proteins, i.e., with similar conductances and life times (some 
early publications: Q [3 |U [T0J [TT] [12)). In their appear- 
ance, they are practically indistinguishable from recordings of 
protein-containing membranes. Although this is surprisingly 
little known, it constitutes a severe problem for the investi- 
gation of biological membranes and protein channels. These 
events will be called lipid channels throughout this manuscript. 
The finding of an enhanced lipid membrane permeability close 
to transitions is particularly important since biomembranes 
seem to exist in a state slightly (about 10-15°) above melting 
transitions. The membrane compositions of fish and bacteria 
adapt upon changes in body temperature, pressure or solvent 
conditions such that this membrane state is maintained (re- 
viewed in lfT3l ). It seems therefore plausible to assuming a 
functional purpose for preserving a particular physical state of 
the membrane. 

The most common ways to measure conductances of mem- 
branes are the black lipid membrane (BLM) technique by 
Montal & Mueller [14] and the patch-clamp technique pio- 
neered by Neher & Sackmann ifTBI . In the Montal-Mueller 
technique an artificial membrane is formed over a small teflon 
hole (normally with a diameter of 50-250 /im). The patch- 
clamp technique relies on the measurement of currents through 
small membrane patches defined by a glass capillary of a 
diameters ranging from 1-10 /im. While the Montal-Mueller 
technique is suited for measuring artificial systems of defined 
composition, e.g., one single protein species reconstituted into 
a single lipid bilayer, the patch-clamp method is mainly used 
for recording currents through biological membranes with 
complex and often unknown composition. Patch pipettes have 
also been used for measurements on synthetic BLMs |[T6l[T7ll . 
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In both techniques, the observed membrane area is much larger 
than that of a protein. The diameter of an ion channel protein, 
e.g. the potassium channel, is about 5nm. This implies that the 
smallest patch in a patch-clamp experiment is already 40000 
times larger than the cross-sectional area of the protein. The 
hole in an BLM experiment is rather 400 million time larger 
than a typical protein channel area (assuming an aperture of 
100 /im). However, the electrophysiological experiment in 
itself cannot tell where along which path the ion currents flow. 
Thus, if one wants to attribute a particular current event to a 
protein one has to assume that at least 99.998% of the sur- 
rounding membrane (including the other membrane proteins) 
are inert in the electrophysiological experiment. Therefore, it 
can generally not be excluded on the basis of the experiment 
itself that during an electrophysiological experiments one also 
finds currents through the lipid membrane. In the Hodgkin- 
Huxley paper [1|, the possibility of leak currents is included. 
If these are structureless small amplitude currents in the back- 
ground this does not represent a major problem. However, leak 
currents represent a significant conceptual problem to protein 
channel recordings if they display a similar signature as the 
protein events. 

The purpose of this review is to demonstrate that in fact 
the ion currents through lipid membranes resemble those 
through proteins to a degree that they become indistinguish- 
able. Their occurrence depends on temperature, lateral tension, 
membrane-associated drugs as anesthetics and neurotransmit- 
ters, pH, calcium concentration, voltage and other thermody- 
namic variables. Here, we review the literature on membrane 
permeability and lipid ion channels, and introduce into the 
thermodynamics of their creation. In this context we discuss a 
possible role for proteins as membrane perturbations that alter 
the state of the lipid membrane. 

2 Macroscopic changes of permeability 
in transitions 

2.1 Experiments 

During the 1960s and 70s, many researchers were investigating 
the permeability of biomembranes for ions. Proteins were 
assumed being the major players. However, accurate controls 
of the permeability of the lipid membrane for ions was needed 
as a control. Hauser et al. [18] studied the permeability of 
lipid extracts from egg yolk or ox brain and found the perme- 
ability for ions being small compared to the permeability of 
biological membranes. The authors concluded that the role of 
the lipid bilayer for permeability is negligible for the case of a 
biological membrane suggesting that the predominant part of 
permeability has to be attributed to proteins. 

However, it was soon shown that the permeability can be 
orders of magnitude larger close to the chain melting regime 
of lipid membranes. Papahadjopoulos et al. were the first 
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Figure 1: Left: Permeability of DPPG vesicles for ra- 
dioactive sodium, 22 Na + , close to the melting transition, 
adapted from ^21/. Right: Permeability of DMPC mem- 
branes for Co 2 +. Courtesy O.G. Mouritsen, cf. /|5]?. The 
dotted profiles are a gide to the eye. The dotted vertical 
line marks the melting transition temperature, T m . 



to demonstrate that the permeability for sodium ions (they 
used radiolabeled 22 Na + ions) increased by at least a factor of 
100 in the phase transitions of dipalmitoyl phosphatidylglyc- 
erol (DPPG) and dipalmitoyl phosphatidylcholine (DMPC) 
in agreement with the phase transitions of these lipids as 
measured by fluorescence changes of embedded markers. The 
permeability curve for DPPG membranes is shown in Fig. [T] 
(left). The permeation time scales in this paper seem extremely 
small (of the order of hours) which might be related to the fact 
that none of the data points was recorded at the phase transition 
directly. The permeation profile for DPPC was found to be 
similar similar (not shown). It was demonstrated that choles- 
terol both abolishes the permeability maximum and the chain 
melting discontinuity. Along the same lines Sabra et al. 
found that the permeability of dimyristoyl phophatidylcholine 
(DMPC) membranes for Co 2+ was drastically enhanced in 
the phase transition regime (Fig. [T] right). These authors also 
demonstrated that the insecticide lindane changes the perme- 
ability. This phenomenon will be discussed in more detail in 
section ROI 

Jansen and collaborators [6| showed that membranes in their 
transition are much more permeable to water (Fig. [2] left). 
Vesicles filled with D2O display a contrast with respect to an 
H2O background leading to enhanced light scattering in an op- 
tical experiment. Permeation for D2O was monitored in a rapid 
mixing stopped-flow experiment. It leads to a mixing of H 2 
and D2O and a loss of scattering contrast. The authors found 
that the permeability of a DMPC membrane changes strongly 
in the phase transition regime such that exchange of water from 
a vesicle is complete after 2ms, the fastest time that could be 
recorded in this experiment. The same finding was reported 
for other lipids as DPPC and distearoyl phosphatidylcholine 
(DSPC), a series of phosphatidyl ethanolamines, phosphatidyl 
serines and phosphatidylglycerols, rendering this study one of 
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the most complete works on membrane permeability. Since 
water permeability typically is discussed in the context of aqua- 
porines |[T9l it is important to point out the the membrane itself 
is highly permeable to water once one is close to a transition, 
which seems to be the case for many biological membranes 
ll20l IT3l . Blicher and coworkers [2| found in fluorescence 
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Figure 2: Left: Permeation time scale for water. Light 
scattering of DMPC vesicles filled with D2O in a H2O 
environment is monitored after a fixed time in a stopped- 
flow experiment. The decrease of the scattering ampli- 
tude close to T m =23.7°C is due to fast mixing of D2O 
and H2O. Adapted from /|6^. The dotted line is a guide to 
the eye. Right: The permeation rate ofRhodamine 6G out 
of vesicles from a mixture of DPPC and DPPG (95:5) in 
the presence of 200mM NaCl. The solid line is the exper- 
imental heat capacity profile for vesicles of the identical 
preparation. Adapted from ^j. 



correlation spectroscopy measurements that vesicles mixtures 
of DPPC and DPPG became significantly more permeable 
for the positively charged fluorescence marker rhodamine 6G 
(R6G) in the phase transition regime (Fig. [2j right). While 
the permeation time scale in the solid lipid phase was of the 
order of many hours, it was 100 seconds and faster in the 
transition regime (the fastest time scale that could be recorded 
in the FCS experiment). The characteristic dimension of this 
marker is approximately lnm. A control experiment with a 
fluorescent-labeled sugar that are larger than R6G did not show 
any measurable permeation on the timescale of hours. This 
demonstrates that 

1. vesicles stay intact and do not rupture in the transition. 
One may confidently assume that this is also the case for 
the experiments shown in Fig. [T] and [2] (left). 2. 

2. sizes of pores responsible for the permeation process 
must be of nm-size. 

The conduction for ions and small molecules described in this 
section was named 'macroscopic' permeability because the ex- 
periments were made on ensembles of vesicles. This kind of 
measurement describes the relation of the lipid melting events 
and the permeability correctly but does not make any statement 



on the molecular nature of any permeation process. The latter 
will be discussed in more detail below in section|U 



3 Techniques for measurement of lipid 
channel activity 

The electrical properties of Lipid bilayers and biological mem- 
branes are measured through small membrane patches defined 
by a hole in a hydrophobic substrate 1 14 1 (black lipid mem- 
branes), the tip of a glass pipette [21 1 (patch-clamp), or more 
recently by a hole in a glass surface on a chip [ 22 1 in a technique 
called 'planar patch-clamp'. Both Montal-Mueller and patch- 
clamp technique have been used to investigate permeability and 
channel formation in pure lipid membranes and are briefly de- 
scribed below. 

3.1 Black lipid membranes / Montal-Mueller 
technique 

In 1972, Montal and Mueller [14] described a method to form 
lipid bilayers formed over a hole in a thin hydrophobic film 
(e.g., teflon) from monolayers. Monolayers are formed on an 
electrolyte surface from a solvent solution, typically containing 
decane, octane, hexane, pentane, chloroform or other apolar 
solvents. The monolayer trough is separated be a teflon film 
into two compartments. A tiny hole in the film (several 10 
fim) created by a needle or by an electric spark slightly is 
located slightly above the monolayer surface. The level of the 
aqueous buffer is now raised and the monolayers fold over 
the aperture in the teflon film as shown schematically in Fig. 
[3] (left). The monolayers merge into a bilayer at the edge of 
the teflon aperture. At these edges one expects considerable 
mechanical tension of quite unclear physics. In order to reduce 
this tension the teflon aperture is normally pre-painted with 
hexadecane, a long-chain alkane that is practically water insol- 
uble. Hexadecane must be considered as a impurification that 
accumulates at the edge of the aperture. It is quite difficult to 
estimate how much hexadecane is dissolved in the membrane. 
For this reason, the composition of the membrane is always 
somewhat unclear even when the composition of the monofilm 
is exactly known. Currents are measured with electrodes in the 
two compartments in a voltage clamp experiment. The capac- 
itance can be measured by continuously changinging voltage 
such that dV/dt = const. The capacitance C m is obtained 
from the capacitive current I c = C m dV/dt necessary to 
charge the capacitor (assuming that the capacitance itself is a 
constant, dC m /dt = 0, which may not necessarily be correct if 
voltage changes the state of the membrane). The measurement 
of the capacitance is required to ensure that only one single 
membrane has formed. Typical capacitances of membranes are 
between 0.6 and 1 uFlcm 2 . 

Membranes forming over a hole appear black in visible light 
due to destructive interference of light reflected from the 
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Figure 3: Two methods to study currents through artifi- 
cial membranes. Left: In the Montal-Mueller technique 
M4V two monolayers merge into a bilayer over an aper- 
ture in a teflon film . Typical aperture size is 50-250 pm. 
Right: One can also generate a membrane patch on a 
glass pipette tip dipped into the surface of a Langmuir 
monolayer trough M6\ 1771/ . Typical pipette tip sizes are 
l-10pm. The drawings were adapted from H23V . The two 
microscopic images on the bottom show the holes in a 
teflon film [24] and a pipette tip used for experiments in 
our lab ( courtesy Katrine R. Laub, NBI). 



two bilayer surfaces. They are called black lipid membranes 
(BLMs) because their formation can be seen by visual in- 
spection under a microscope. The typical hole diameter in 
a Montal-Mueller experiment is 50 to 250 /im (Fig. [3] left 
bottom). The disadvantage of this technique is the presence of 
solvent impurities and the large tension at the edge of the teflon 
aperture. Furthermore, due to the large diameter of the hole the 
membranes are often quite unstable. However, this technique 
is often used to study lipid membrane permeability but also the 
conductance through channel proteins and peptide pores. 

3.2 Patch-clamp technique 

The patch-clamp technique has been introduced by Neher and 
Sakmann 1 15 1 and was described in detail by Hamill et al. [21 1. 
It was first used to study channel proteins in biological cells (in 
particular the acetylcholine receptor in frog muscle cells). A 
small pipette is placed on a cell surface while applying slight 
suction. Currents and capacitances are measured by using two 
electrodes, one being located in the pipette and the other one 
in the cell. The patch pipettes have much smaller diameters 
of order 1-10 pro. (see Fig. [3] right bottom). The obvious 
advantage over the Montal-Mueller technique is that one can 
investigate intact cells and on does not have to use organic 



solvents. The interface of the patch pipette and the suction is a 
source for perturbation GBI . 

Patch pipettes can also be used to study BLMs lfl6llT7l . The 
tip of a glass pipette filled with electrolyte solution is placed on 
a buffer surface. Lipids dissolved in an organic volatile sol- 
vent (we use mixtures of hexane and ethanol) are brought into 
contact with the outer surface of the patch pipette. When the so- 
lution runs down the pipette a membrane spontaneously forms 
at the pipette tip as schematically shown in Fig. [3] (right). Af- 
ter membrane formation as evident in the capacitance measure- 
ment slight suction is applied. In this technique no long chain 
alkanes are used and one can assume that after some equilibra- 
tion time the membrane at the pipette tip is practically solvent 
free (W. Hanke, personal communication). Hanke and cowork- 
ers ifTTIl also describe a variation of this technique called tip- 
dipping where the membrane patch is generated from dipping 
the pipette tip into monolayer surfaces. 

4 Lipid channels 

The evidence for pores (channels) within membrane proteins 
selectively passing ions rests partially on the observation that 
the currents through membranes are quantized (Fig|4^). The 
fact that the currents are small and correspond to aqueous 
channels of the dimensions of a few A suggests that the con- 
duction of ions occurs along objects of molecular dimension. 
They have been attributed to particular proteins called ion 
channels. The literature about these proteins is huge (e.g., 
[26 1). It seems therefore surprising that one can find similar 
quantized current events in pure lipid membranes that do not 
contain proteins. Fig. [4^ shows the historical experiment of 
Neher & Sackmann [15| on frog muscle cells containing the 
acetylcholine receptor in the presence of an agonist molecule. 
Fig. [4J5 shows the currents through a lipid membrane from 
a mixture of dioleoyl phosphatidylcholine (DOPC) and di- 
palmitoyl phosphatidylcholine (DPPC) having a chain melting 
transition close to the experimental temperature. The melting 
temperature is important since it has been shown in section [2] 
that membranes are most permeable in their melting transition. 
It is obvious that the two traces display very similar features 
both regarding amplitude and time scales of the currents. In 
the following we will characterize the lipid channels in more 
detail. 



4.1 Early experiments 

To the best of our knowledge the first channel events recorded 
in pure lipid membranes are from Yafuso et al. in 1974 [7], two 
years before the famous paper by Neher and Sakmann on the 
acetylcholine receptor shown in Fig. [4^. Their experiments on 
bilayers from oxidized cholesterol are shown in Fig. [5] (left). 
They found spontaneous conductance changes and multilevel 
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(a) acetyl choline receptor in frog muscle cells 
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Figure 4: Comparison of conductance events through bi- 
ological membranes containing channel proteins and a 
pure lipid membrane in its melting regime, a. Record- 
ing of the acetylcholine receptor in frog muscle cells from 
the classical paper of Neher & Sakmann i T/51/ (Ringer 
solution in the pipette). b. BLM experiment on a 
DOPCDPPC membrane in 150mM KCl, 60mV in its 
transition regime at 19° C. Figure adapted from [27]. 



conductance states in an unmodified synthetic lipid film (of 
somewhat exotic nature). While the lipid system used does 
not allow to draw conclusion with respect to biological mem- 
branes, the experiments quite clearly are a proof-of-principle 
in so-far as they demonstrate that lipid films alone can generate 
phenomena that are typically attribute to proteins, peptide 
pores or other additives. The original recording was made 
with an analogue chart recorder in which the pen position was 
changed by rotation around an axis. This has been digitally 
corrected by us (see the original paper for a reference to the 
raw data). A can be seen in Fig. [5] (top) a change of the 
voltage leads to an overall larger membrane current the shows 
quantized steps. 

Another important early paper showing this phenomenon 
is authored by Antonov and collaborators in 1980 |8|. They 
showed that one obtains quantized currents DSPC membranes 
one obtains quantized currents of in the pA regime at 100 
mV. As in the paper by Yafuso et al. Q, they found that 
current traces may display several steps as shown in the current 
histogram in Fig. [5] (bottom). In their paper they state: "We 
suggest that these channels could conduct the transmembrane 
ionic current in biological membranes without the involvement 
of peptides and proteins", thus making the first statement 
indicating that this could be a general phenomenon in biomem- 
branes. 

Other authors have found the same phenomenon. Fig. [6] (top) 
shows results from Yoshikawa et al. |23l . They described 
quantized current fluctuations through DOPC membranes con- 
taining some cholesterol that are of order 1 pA (corresponding 
to a conductance of 50 pS). The authors found similar results 
for membranes in BLMs, on patch pipettes and on membranes 
deposited on filter paper. As other authors, they also discuss 
the relevance of this finding for the interpretation of currents in 
biomembranes and suggest that currents attributed to proteins 
can also arise from lipid membrane fluctuations. 

Fig. [6] (bottom) shows quantized membrane currents in a 
membrane made from POPE:POPC=7:3 mixture (from |28|). 
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Figure 5: Top: First channel recordings in a lipid mem- 
brane made of oxidized cholesterol [7]. The bottom line 
indicates the voltage. The dashed lines indicate the volt- 
age changes that give rise to large capacitive currents in 
the current recording. For constant voltage current fluc- 
tuations of the order of 10-20 pA occur. The data were 
digitally adapted from the original Figure and numeri- 
cally corrected for the distortions caused by the analogue 
chart-recorder. Bottom: Quantized currents through syn- 
thetic DSPC membranes close to their transition in 1M 
KCl at WOmV W- 



30 min 



Figure 6: Top: Current fluctuations in a mixture of DOPC 
and cholesterol (about 9:1) at 20mV in 1 M KCl A23\I . 
Bottom: Conductance fluctuations of a membrane made 
of a POPE:POPC=7:3 mixture in 150mM KCl l[28\l . 



The currents display a similar conductance as in the above 
experiment. The authors discuss the similarity of the lipid 
conduction steps with those attributed to proteins channels and 
state: "This observation demonstrates the need for caution in 
interpreting conductance changes, which occur following ejec- 
tion of channel- containing vesicles near a membrane." This 
refers to the standard reconstitution procedure of generating 
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Figure 7: Currents through a purely synthetic diphytanoyl 
phosphatidylcholine membrane ( 1M KCl unbuffered, pH 
ca. 6.5). Left: Currents recorded at 77mV in differ- 
ent scalings. Right: Current histogram showing several 
steps. From Kaufmann et al. ISSS - 



synthetic membranes containing proteins. In particular, the 
authors state that currents attributed to proteins could in fact be 
due to lipid pores. 

Kaufmann et al. [29 1 showed further interesting results on 
channel appearance in purely synthetic membranes. Their pa- 
per appeared as a booklet (available from the website given in 
the reference list and from the University Library in Gottingen, 
Germany) and gives considerable insight into the physics of 
the pore formation process that goes along the line of argument 
in section 5.2 A nice example for quantized currents through 
membranes is given in Fig. |8](top) 0. At least 4 steps in the 
current recording through a DOPC:DPPC= 2:1 mixture can be 



seen. Each of them corresponds to a conductance of about 300 
pS. The experimental temperature of 19°C is close to the tran- 
sition maximum of this lipid mixture. Fig. [8] (bottom) shows 
a long trace from a BLM made of a DCi 5 PC:DOPC=95:5 
mixture in the phase transition regime (31.5°C, lOOOmV) |30|. 
The above experiments show that there is quite a range of pos- 
sible pore sizes and conductances found for lipid membranes. 
We focus here on literature demonstrating quantized current 
events through lipid membranes. There is additionally a com- 
plete literature on reversible electrical breakdowns of artificial 
membrane conductance at higher voltages [31, 32, 33], which 
are connected to electroporation. We believe, that lipid channel 
formation and electroporation are related phenomena. 

4.2 Dependence of channel formation of the 
melting transition 

As mentioned in section [2] the overall membrane conductivity 
is maximum in the transition regime. Consistent with this find- 
ing the likelihood of finding lipid channel events is maximum 
in the lipid melting regime. Below we show three experiments 
that demonstrate this. Boheim et al. [9] reported quantized cur- 
rent events in synthetic 1,3-SMPC membrane at the phase tran- 
sition temperature of about 29°C (Fig. [9j. The discrete conduc- 
tance steps are not well visible in the figure, but were obvious 
in the original data trace (Wolfgang Hanke, personal commu- 
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Figure 8: Top: Current fluctuations in a 
DOPC:DPPC=2:l mixture at 19° C (150mM KCl, 
60mV). One can resolve at least 4 quantized steps with a 
current amplitude of about 20 pA corresponding to about 
to a conductance of about 300pS. From Bottom: 
Current fluctuations in a DC\§PC:DOPC=95:5 mixture 
at 31.5° C (lOOOmV), which is the transition regime of 
this mixture. Quantized currents with an amplitude of 
15pA can be seen. From [30]. 



nication). They write: "Apparently, the transition of a bilayer 
membrane from the liquid-crystalline to the solid state and vice 
versa is associated with spontaneous fluctuations in membrane 
conductance." Antonov et al. demonstrated for both DSPC [8 1 
and DPPC 1 34 1 that they found quantized currents in the phase 
transition regime while the membrane seemed inert outside of 
this regime. Fig 10 (left) shows DPPC membranes above the 



transition at 50°C (top), in the melting regime at 43°C and be- 
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Figure 9: Currents through a 1,3-SMPC membrane with 
a transition at 29° recording during a temperature scan. 
At the transition temperature quantized current events 
were found (not well resolved in the figure but well de- 
scribed in the text). From [9 J. 
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Figure 10: Currents through DPPC membranes with T m 
around 41° C. Left: heat capacity profile of DPPC. Right: 
Currents through DPPC membranes above (50° C), close 
to (43° C) and below (35°C) the transition. Only close 
to the transition quantized currents with discrete steps in 
the current histogram were seen. From [34]. 



low the transition at 35°C. The current fluctuations in the tran- 
sition are quantized and are of unusual magnitude (nA regime 
indicating very large pores). 

Transitions in single lipid membranes occur over a quite nar- 
row temperature interval of less than one degree which makes 
it difficult to adjust the experimental temperature accurately. 
The heat capacity displays a very large amplitude at maximum. 
As shown below (section |5.5[ ) that the permeability is closely 
related to the heat capacity, and that fluctuations in area are 
very large. Thus, in transitions of single lipid membranes they 
display a tendency to show very large pores and a pronounced 
tendency to rupture. Those black lipid membranes are inher- 
ently unstable and break easily. 

For the above reason, other authors prefer to work with lipid 
mixtures. They display broader melting profiles with smaller 
fluctuations at maximum. These membranes are more stable 
and due to their wider melting regime not as temperature sen- 
sitive. This approach has been taken by [24]. They used a 
DOPC:DPPC=3:l mixture (150mM KC1, 40mV) that displays 
a broad melting profile with a maximum around 17°C.At the 
heat capacity maximum many discrete current steps of about 
30pA were found (Fig. 
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, corresponding to a conductance of 
i=3 750 pS. Above the melting temperature at 30 °C (in the fluid 
phase regime) no discrete current steps were seen. 

4.3 Specificity of lipid channels 

In their seminal paper from 1980, Antonov et al. |8 1 compared 
quantized currents through membranes made of the zwitteri- 
onic DSPC with those through positively charged membranes 
of a DSPC analogue in which the P-0 in the head-group 
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Figure 11: Currents through BLMs of a 
DOPC:DPPC=3:l membrane (from /|2]?). Top: heat 
capacity profile with a maximum at 17° C. The two 
dotted lines represent the temperatures at which the 
currents in the bottom panel were recorded. Bottom, left: 
Current steps could only be observed at the transition 
temperature {17° C, 150mM KCl, 40mV) but not in the 
fluid phase regime at 30° C. Bottom, right: Current 
histograms. A baseline current was subtracted. 



was replaced by a P-CH3 rendering the membrane positively 
charged. They found that the cationic membranes were anion 
selective while zwitterionic membranes did not display any 
selectivity for charges (unfortunately, no experimental details 
were given). 

In [34- 1, the same group of authors compared the conductances 
of membrane pores for different cations. They recorded the sin- 
gle channel currents through DPPC membrane in the presence 
of LiCl, NaCl, KCl, RbCl and CsCl. In a voltage regime from 
-80 to +80 mV they found a linear current-voltage relationship 
from which a conductance gi can be calculated (Ii = gi ■ U, 
where i is the index for the ion species). The details of the 
experiment are given in Fig. 



12 The authors found that the 



value of the conductance decreases in the sequence: 

9Li > 9Na > 9K > 9Rb > 9Cs 

Antonov and collaborators discussed this sequence in con- 
nection with kosmotropes (water structure enhancers) and 
chaotropes (water structure breakers) that are closely related to 
the Hofmeister series. This series involves a ranking of ions 
in respect to their effect on various materials in electrolytes 
(e.g., proteins) and is related to the interaction of the ions with 
water. The Hofmeister series has been discussed in connec- 
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Figure 12: Currents through pores in DPPC membranes 
close to the chain melting temperature (43°). Top: Cur- 
rent traces in the presence of LiCl at different voltages 
lead to a linear current-voltage relationship. Bottom: 
Current-voltage relationships in the presence of LiCl, 
NaCI, KCl, RbCI and CsCl yield different conductances. 



tion with peptide pore conductances (35). According to 
the sequence from chaotropic to kosmotropic for the above 
ions is Li + > Cs + > Na + > K + while it was given as 
Na+ > K+ > Rb+ > CS+ by M. ChaplirQ consistent with 
the conductances in the above experiment. It is very likely 
that such correlations exists, rendering the lipid membrane 
moderately ion selectiv^] 



5 Theoretical considerations 
5.1 Pore models 

The remarkable finding in electrophysiological lipid membrane 
experiments is that ion permeation expresses itself as quantized 
current events. This all-or-nothing feature is unexpected since 



1 See: wwwl.lsbu.ac.uk/water/hofmeist.html 

Selectivity is considered here to be a conduction preference expressed in 
different permeation rates rather than an exclusive all-or-nothing conduction 
phenomenon. 



one might expect the membrane defects form through stochas- 
tic processes that favor broad distributions of pore sizes. The 
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Figure 13: Two representations for transport across a 
membrane. Top: Single file transport. Water molecules 
align in a pore. From Finkelstein 4571/ . Center: 
Schematic drawing of membrane pores: a hydrophobic 
pore. Bottom: a hydrophilic pore. From Glaser et al. 

ess. 



permeability of membranes has been imagined to occur via 
three mechanisms: 

1 . diffusion of solutes through the hydrophobic core of the 
membrane 

2. transport along water files consisting of single H 2 
molecules (FigfT3]left) 



3. hydrophobic pore formation (Fig 1 3 right top) 



4. hydrophilic pore formation (Fig 13 right bottom) 



Mechanism 1 was proposed for oil-soluble substances like 
anesthetics already in the 19 th century 139) . However, for 
ions or polar molecules as water this is an unlikely mechanism 
considering the very low solubility of such substance in the 
hydrophobic core membrane core. The free energy of a charge 
in a dielectric medium is proportional to the electrostatic po- 
tential W = q/eoer where q is the charge, en. is the vacuum 
permittivity, e is the dielectric constant, and r the distance 
from the charge. The dielectric constant in water is about 80 
while it is 2-4 in the membrane interior. Thus, the electrostatic 
free energy is different by about a factor of 20-40 between 
membranes and water. Due to Debye-Hiickel ion screening 
this effect is even larger in the aqueous electrolyte solutions. 
Transport mechanism 2 has been described by Levitt ATI and 
Finkelstein ll37ll . see Fig. [13] (top). Finkelstein defines: "By 
single-file transport of water one means that water molecules 
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Figure 14: Molecular dynamics simulation of a POPC 
membrane with a tansmembrane voltage of2V. Top: Hy- 
drophobic prepare in an early stage of the simulation 
(less than Ins). Bottom: Pore with dimensions of the or- 
der of 1 nm after about 50ns. Adapted from Bockmann et 

al. m. 



within the pore cannot pass or overtake each other." In this 
context the term "no-pass transport" has been used. Ions would 
be part of the file and can be transported in sequence with the 
water molecules. The single file mechanism might explain the 
quantized nature of currents because there is a minimum size 
for the water file defined by the dimensions of H2O and the 
ions. 

Pore formation is the most popular mechanism proposed 
for membrane transport. Glaser et al. |[38l suggested both 

center) and hydrophilic 



hydrophobic (mechanism 3, Fig 13 



pores (mechanism 4, Fig 13 bottom). The hydrophobic pore is 
probably related to the formation of a water file as described 
in (2) and does not involve any rearrangement of lipids. The 
hydrophilic pore involves a change in the geometry of the lipid 
head groups. The elastic free energy of such a pore can be 
calculated assuming typical values for the elastic constants. 
Stable pores with a diameter of ~ 1 nm were suggested. On 
theoretical grounds it was further proposed that transmem- 
brane voltage can stabilize large pores at high voltage |42|. 
This might be important in connection with the phenomenon 
of electroporation ll33l . 

Recent molecular dynamics simulations of membranes ex- 
posed to large voltages (~ 2V) suggest that both single file 
and hydrophilic pore formation can occur. Fig. 14 shows a 



POPC membrane at two stages during the simulation after a 
voltage of ~ 2V was applied. In Fig. [14] (top) a single water 
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Figure 15: Schematic drawing of the chain of events in 
the molecular dynamics simulation of POPC membranes 
shown in Fig. 14 After the application of a voltage of 
2V across the membrane, first a "prepare" opened that 
is both consistent with the water file shown in Fig. [75] 
(top) and the hydrophobic pore in Fig. ^13^(centerj. Sub- 
sequently, a hydrophilic pore with a radius of about 0.5 
nm formed. From Bockmann et al. [40]. 



file formed along which ions could flow. This was called a 
'prepore' [40] and occurred after less than 1 ns of simulation 
time. After ~ 50ns a large pore developed long this defect that 
had a diameter about 1 nm (Fig. [14] bottom). The chain of 
events in the simulation is shown in Fig. 
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First, a water file 
This is consistent 
(Fig. 13 center). 



(prepore) forms along which ions can flow, 
with the hydrophobic pore suggested by [32 
Subsequently, the lipids rearrange and form a hydrophilic pore 
(Fig. [T3] bottom). 

If we assume that this is the likely chain of events then we 
would conclude that the single file (prepore) defines a minimum 
size of a conducting event in a membrane. The hydrophilic pore 
seems to display a stable radius of the order of 0.5 nm. This is 
consistent with the magnitude of the currents measured in the 
patch-clamp and the Montal-Mueller techniques. Remarkably, 
this is also the typical size of the aqueous pores proposed for 
channel proteins. 

Assuming a pore of radius r, Glaser et al. 11381 calculated a free 
energy of pore formation of 



AF(r) = 2ir-ir + irr 2 (a + aU 2 ) 



(1) 



where 7 is the edge tension of the pore, a the lateral tension 
and U the voltage, a is a coefficient related to the capacitance 
of the membrane. The free energy therefore consists of three 
contributions: the energy of the pore interface, the work that 
has to be performed against an external tension a that is pro- 
portional to the pore area, and a correction term containing 
voltage that is also proportional to the pore area. One should 
assume that the probability of pore formation is related to 
exp(— AF(r)/kT) and the conductance of one pore is in a 
simple manner related to the pore radius r. 

The above equation does not yet contain a particular the- 
ory for the effect of the phase transition, which is obviously of 
major importance as judged by the experimental data shown in 
sections [2]and|4] This will be provided below. 
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5.2 



Entropy, 
forces 



thermodynamic variables and 



In 1910, Albert Einstein wrote a remarkable paper on critical 
opalescence B3l in which he considered the phenomenon of 
maximum light scattering of liquid mixing close to a critical 
point (called "critical opalescence"). This paper contains a 
deep reflection of the nature of the second law of thermodynam- 
ics and is at the basis of all non-equilibrium thermodynamics 
and fluctuation theory. The second law basically states that 
the most likely state of a system is most probable. This were 
without content if one could not extract more information from 
this law without making assumptions on the molecular detail 
of a system. Einstein showed in fact that one can make general 
statements about the nature of fluctuations. His theory for 
critical opalescence is similar to the thermodynamic theory of 
the pore formation process outlined below. The application of 
this theory for pore formation was already proposed by Nagle 
& Scott it) and Kaufmann et al. 11291 . 

Let us consider the entropy as a function of the extensive 
thermodynamic variables a (e.g., energy E, volume V, area A, 
charge q, ...) and define the deviations from the equilibrium 
value £,i = ai — a* o we can expand it into a Taylor series: 



thermal 
equilibrium 




variable 



Figure 16: The entropy as a harmonic potential and the 
corresponding fluctuations of the variable £ that could be 
internal energy, volume, area etc. The slope of the poten- 
tial at £ is the thermodynamic force, while the curvature 
of the potential at £o is proportional to the corresponding 
susceptibility (cf. eqs. [7 10 and 12 1. 



is given by 07) 



S = S +J2 



as 



6 



d 2 S 



(2) 



k 



The fact that the entropy is at maximum in an equilibrated sys- 
tem implies that the linear terms (dS/d£i)o are zero and thus 
we can approximate the entropy for small fluctuations with 



where (...) denotes the statistical mean, which implies 



/ d 2 S 



_k_ 



(6) 



(7) 



where = —0.5(d 2 S/d£i£j)o are positive constants 



The expansion of the entropy into a quadratic function of 
the intensive thermodynamic variables is at the basis of the 
complete discipline of linear non-equilibrium thermodynam- 
ics as developed by Onsager I44ll45l and Prigogine [46 1 . The 
probability to find a particular pair of fluctuations, £j and £j, is 
given by 



= Sq — 9ij€i€j (3) This corresponds to the mean square amplitude of the fluctu- 
ij ations or the width of the Gaussian distribution. This implies 

that the half width of the distribution is related to the curvature 
of the entropy potential, a fact that we will make use of below 
in order to derive the relaxation times of fluctuations. 



The thermodynamic forces are defined as the derivatives of 
the entropy with respect to the extensive variable. 



Xi 



dS 



9V$> 



(8) 



P&,£i)=Poexp 



5(6, fr) 
k 



P exp 



or for one single fluctuation it is given by 



(4) 



(5) 



This is completely in analogy to the definition of forces in 
classical mechanics (as the derivative of a potential) but more 
general since it includes extensive variable as the energy and 
charge, but also the particle numbers . Further, the potential 
under consideration here is the entropy. 
Some of the thermodynamics forces are 
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with the conjugated variable E (interal energy) 



which corresponds to a Gaussian distribution of states with 
Po = y/gu/n k. T 

The thermodynamic average of the product of the fluctuations — with the conjugated variable V (volume) 



T 
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n 

T 
f 

ti 

T 



with the conjugated variable A (area) (9) 

with the conjugated variable q (charge) 

with the conjugated variable rij (particle number) 



where T is the temperature, p is the pressure, II is the lateral 
pressure, \& is the electrostatic potential, and fXj is the chemical 
potential of species j P7ll46l . E, V, A, q and n^ correspond to 
the fluctuating variables £j. 
The above implies that 



k 9l 



(10) 



meaning that the fluctuations are related both to the curvature 
of the entropy potential of the change of the forces upon the 
change of a variable. Eq. [I0]has been strictly derived in 11471 . 

5.3 Susceptibilities 

The susceptibilities or response functions are functions like the 
heat capacity, the compressibility or the capacitance. They are 
given by the derivative of an extensive variable with respect to 
an intensive variable. E.g., the heat capacity at constant volume 
is given by 

Now, = E — (E) is one of the quantities that can fluctuate, 
while dT = —T 2 d(l/T) is related to the conjugated thermo- 
dynamic force, Xe = l/T. In particular, 



T*\dX E ) v T*\d£ E/ , 
Another way of writing this is 

(E 2 )-(E) 2 



Cy = 



kT 2 



kT 2 
(12) 



(13) 



and the isothermal area compressibility 

_ (A 2 )-(AY 
A kT 



1 

A 



dA 



(16) 



Eqs. 13 to 16 are linked to the fluctuation theory (see, e.g.. 



|47|). The principal underlying idea is that small forces gener- 
ated thermal motion generate local fluctuations in the extensive 
variables. The coupling of the susceptibilities to the fluctua- 
tions is important and will play an crucial role below. 



5.4 



Coupling of compressibilities with the heat 
capacity 



that for transitions in membranes 



It has been argued by 
one finds that 

AV(T) = lv ■ AH(T) (17) 

where AV(T) is the temperature dependent volume change in 
the transition regime, and AH(T) is the corresponding change 
in enthalpy. This means that volume and enthalpy change in a 
proportional manner with a proportional constant ^y = 7.8 • 
10~ 10 m 3 /J for many lipid systems including biological mem- 
branes ll49l . More indirectly, it has been argued that a similar 
relation is also true for area changes in melting transitions of 
membranes: 

AA{T) = 1A ■ AH{T) (18) 

with 7^4 = 0.89m 2 /J ll48l . The two above relation together 
with eqs. 15 and 16 imply that the fluctuations of enthalpy are 
proportional to the fluctuations in volume or area (i.e., enthalpy, 
volume and area are no independent variables), and we obtain: 



a v _ lyT a 

LXKr^p — — — - — L\Cp 



and 



A 



(19) 



This means that both volume and area compressibility are pro- 
portional to the heat capacity changes in membrane transitions. 
In particular it implies that membranes are very soft in their 
transitions. This change in compressibility can easily of order 
10-fold at the transition maximum. This has profound conse- 
quences for the pore formation process. 



This implies that the heat capacity is proportional to the 

fluctuations in internal energy. This relationship could 5.5 Permeability 

also be derived in a simpler manner by just differentiating 
the statistical mean of the internal energy that is given by 
( E ) - EiE i exp{-E i /kT)/'E i exp{-E i /Kr). Thus, the 
fluctuation relations would be strictly true even if the assump- 
tion of a harmonic entropy potential is not made. Similar 
relations can be written for the heat capacity at constant 
pressure 

_ (H 2 ) (H) 2 



A couple of authors have discussed area fluctuations as the ori- 
gin of membrane pores (or lipid ion channels), e.g. (3] [50] |6] 
I21I5T1 but in particular Nagle & Scott [4] and Kaufmann \29\ 
(much of the line of argument in sections 5.2 and 5.3 was given 



in the latter reference and applied to membrane pores). 
There are in particular two views: 



kT 2 

the isothermal volume compressibility 



1 

V 



dV 
dp 



(v 2 ) (vf 

V kT 



(14) 



(15) 



1 . In the transition one finds phase coexistence and the like- 
lihood of finding a membrane pore is at maximum at the 
domains boundary. Therefore, the permeability is pro- 
portional to the length of the domain interfaces. This 
view has been taken by Papahadjpoulos et al. (3) and 
Mouritsen and coworkers I50ll52ll. 
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2. In order for a membrane pore to be created work has to 
be performed that is proportional to the lateral compress- 
ibility 012 EH. 

These two views are conceptually similar but not identical. 
Clearly, the domains occur in temperature regimes with larger 
area fluctuations. More precisely the fluctuations in domain 
sizes are responsible for both the large heat capacity and the 
large compressibility. Domain sizes in lipid transitions depend 
on the cooperativity of the melting process. In computer sim- 
ulations, cooperativity is usually modeled by an unlike nearest 
neighbor interaction that is effectively a free energy penalty 
for the creation of a domain interface. For this reason, more 
cooperative transitions with narrower and more pronounced 
heat capacity peaks display larger domains at the melting 
temperature (50% of the lipid matrix is in a liquid state) and 
consequently a smaller overall domain interface than less 
cooperative transitions. Thus, large fluctuations imply that 
domains are large and that the total domain interface is smaller. 
The data shown, e.g., in Fig. [2] (right) indicate, however, that 
the permeability is about proportional to the heat capacity. 
This favors the second view. The discrepancy between the 
two views would disappear if one would not consider domain 
interface is one-dimensional lines but rather as narrow areas of 
large fluctuations. However, under these circumstances, one 
would have to take both length and thickness of the interfacial 
region into account. For these reasons we follow the argument 
that the overall fluctuations are responsible for the increase 
in permeability. Domain interfaces are discussed again in 
section l8Tl and it is shown that the fluctuations within different 
domains may also be different. 

Nagle and Scott [4| calculated the free energy (the work) 
necessary for the formation of a pore of area A A in a membrane 
with total area to be 



AF = 1 1 



2 fcrp A. 



(AA) 2 



(20) 



with the lateral compressibility of the membrane. This cor- 
responds to Hooke's law for a mechanical spring. The likeli- 
hood to find a pore of size AA is given by 



P(AA) = P exp 



AF 

W 



(21) 



which implies that the likelihood for a pore is the same when 
(AA) 2 1 Kt is the same. This means one finds pores of twice 
the area if the compressibility is four times larger, or that the 
likely pore size scales with 



(AA) 2 oc 



(22) 



As a next step it was assumed that the permeability P of the 
membrane is related to the size of the defects such that one can 
make a series expansion of the permeability as a function of 
area change AA: 



where Pq is the permeability of the pure phase in the absence of 
fluctuations. In the transition one expects fluctuations AA that 
are both positive and negative and the linear term is neglected, 
yielding 



P = Pn 



C 2 {AA) 2 £ Pf„ 



G*2 fcrp 



(24) 



This equation reflects that one expects larger pore areas AA or 
a larger number of pores if the compressibility is higher. Nagle 
& Scott |4| argued that one obtains the same law if one takes 
into account that one finds distributions of pore sizes. 

In eq. 19 we have derived that for lipid membranes in tran- 



sitions the compressibility changes are proportional to the ex- 
cess heat capacity. Therefore, we finally obtain for the perme- 
ability P: 

P = P + aAcp (25) 

where a is a constant. To summarize, one expects that the mean 
pore size and the pore number is larger when the compressibil- 
ity is higher. This leads to a higher permeability. For the above 
calculation, it does not make a difference if a pore has twice the 
area or if one simultaneously finds two pores. Obviously, in the 
above derivation the free energy of the interface between pore 
wall in water has been neglected. This can only be justified by 
experiment. The data in Fig. |2](right) which yield a satisfactory 
proportionality between heat capacity and permeability suggest 
that in fact the free energy of a pore is dominated by the area 
work as assumed in eqs 24 and 25 This neglects the edge term 



P = P + dAA + C 2 (AA) 2 



(23) 



of the pore which was contained in eq. U\ The effect of voltage, 
that was heuristically included in eq. His also not considered 
here because a satisfying description of the effect of voltage on 
the melting transition is yet missing, even though experimental 
evidence suggests that this effect must be significant (see also 
section [6~4] >. 

6 Influence of changes of variables on 
the lipid channels 

From the above derivations it is clear the changes in all inten- 
sive variables should have a predictable influence on the perme- 
ability when their effect on the melting transition is known. Ac- 
cording to the list given in eq. |9]this includes temperature, hy- 
drostatic pressure, lateral pressure, electrostatic potential (volt- 
age) and the chemical potentials (e.g., protein, proton or cal- 
cium concentrations). In this paragraph it will be shown that 
all of these dependencies exist. This is especially so if the lipid 
membrane is close to a melting transition (or any other kind of 
transition). It was shown by I1201 (findings are also reviewed 
in |[T3l ) that biomembranes are in fact often found close to a 
melting transition in their membranes under physiological con- 
ditions (typically about 10-15° above the heat capacity maxi- 
mum of a relatively broad transition). Therefore, any variable 
that potentially changes the transition will display an influence 
on the occurrence of lipid ion channel events. That this is in 
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fact the case is shown below for a number of thermodynamic 
variables. 

6.1 Temperature 

The previous section justified why the formation of lipid chan- 
nels is most likely in the melting transition of lipid membranes. 
This implies that temperature must be an important variable. 
Let us consider the following very simple mass action scheme 
for a transition at constant pressure 



K 



[F] 
[G] 



exp 



AG 
RT 



exp 



AH 



RT 



TAS\ 



(26) 



where [F] is the concentration (or fraction) of fluid lipids and 
[G] is the concentration of gel lipids, respectively. In the fol- 
lowing we will assume that the state of the membrane and its 
heat capacity depends on AG. For AG > we find a larger 
fraction of gel phase, while for AG < the membrane is 
mostly in its fluid phase. The melting temperature, T m , shall 
be defined as the temperature where [G] = [F] , or K = 1 and 



AH -T m AS = 



AS 



AH 



AG = AH 



T 



(27) 



Thus, obviously a change in temperature will lead to a change 
of the free energy difference between gel and fluid phase. 
Since at the melting point, T m , AG is zero, the fluctuations 
are maximum because the free energy necessary to move the 
system from one state into the other one is zero, and thermal 
noise is sufficient to do so. Therefore, the likelihood to form 
pores and the permeability must be strongly influences by 
temperature. This was shown for the macroscopic permeability 
experiments in secti on [2] |[3l |6] EJ and for the lipid channel 
formation in section |4~2| | 34l l2l . 

The effect of temperature could be expressed as being a 
'thermo-sensitivity' of channel formation in pure lipid mem- 
branes. 



6.2 Pressure and lateral tension 

In the following, we will express the state of the system by the 
Gibbs free energy difference between fluid and gel state. We 
will show that it depends not only on temperature but also on 
all the other intensive thermodynamic variables. 

The enthalpy is a function of pressure: H = E + pV. If 
AH n is the enthalpy difference between fluid and gel phase at 
a pressure of 1 bar, AV the respective volume difference, and 
Ap the pressure difference with respect to 1 bar, we obtain 



T m {Ap) 



AH 
~AS 



AH + ApAV 
AS 




\ . 



SOmS 

Figure 17: Channel events in a membrane made of soy- 
bean phosphatidylcholine in 1M unbuffered NaCl, glass 
pipette tip diameter ca. 5fim, at 500 mV. The top trace (A) 
recorded before application of suction does not show any 
current fluctuations. B. The bottom traces were recorded 
with suction applied. The change in membrane area was 
monitored by recording the membrane capacitance that 
increased from 260 to 510 pF. From [29]. 



AiJ (l + 7 vAp) 
AS 

T m (l+ lv Ap) 



AT m (Ap) = j v T m Ap 



(28) 



(29) 



where 7^ = 7.8 • 10~ 19 m 3 /J is the parameter defined in eq. 17 
Similarly, the shift of the melting transition by lateral pressure 
(or tension) is 



AT m ,(An) = 7A T m An 



(30) 



18 1. It has been shown by 



with = 0.89 m 2 /J (see eq. 
Kaufmann et al. [29] that the application of suction on a 
patch pipette can generate channel events. Suction changes the 
lateral tension/pressure of the lipid membrane and thus has an 

Fig. 



17 shows 



influence on the state of the lipid membrane, 
suction induced channel formation in soybean phosphatidyl- 
choline in 1M NaCl at 500m V. 

The influence of pressure on channel appearance could be 
called a 'mechano-sensitivity' of channel formation in pure 
lipid membranes. 
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Figure 18: Effect of the anesthetic octanol on the melting transition of a D0PC-DPPC=2:1 mixture in 150mM 
KCl at 19° C. Left: Heat capacity profiles in tlie presence of increasing quantities of octanol demonstrate the shift 
of the transition caused by the anesthetic drug. The vertical dashed line indicates the experimental temperature. 
Center: Membrane current recorded at 210 mV demonstrates that octanol both decreases the frequency of current 
events and the amplitude until currents are completely 'blocked'. Right: Current histograms of the center traces. 



6.3 Anesthetics, neurotransmitters and other 
drugs 

Anesthetics can lower melting temperature by AT mXA due to 
a very simple effect called melting point depression [53|. It is 
described by 



AT„ 



AH 



X A 



(31) 



where xa is the molar fraction of anesthetics dissolved in the 
membrane, AH is the melting enthalpy of the membrane, and 
T m is the melting temperature. To arrive at this law one has 
to assume that the anesthetics molecules dissolve ideally in the 
fluid lipid membrane, while they are insoluble in the gel mem- 
brane. The success of eq. 3 1 in describing the effect of anes- 



thetics on melting transition [53] does not only justifies this 
assumption but also explains the famous Meyer-Overton rule 
that states that the effectiveness of an anesthetic is proportional 
to its solubility in oil (the membrane core), independent of its 
chemical composition I54ll55ll56ll . 

Fig. [18] shows the effect of the anesthetic octanol on mem- 
branes made of a DOPC-DPPC mixture. In the absence of oc- 
tanol, the experiment is performed at the heat capacity max- 
imum where many channel events can be found. Increasing 
amounts of octanol shift the transition maximum towards lower 
temperatures and thereby lower the heat capacity at the ex- 
perimental temperature. As a consequence, the frequency and 
amplitude of the current events decreases until channel events 
completely disappear. This effect is also discussed in detail in 
[2,51]. If the membrane is recorded at a temperature below that 
transition, anesthetics can drastically increase the membrane 
permeability because they move it into the transition regime 
0. 

Sabra et al. J5 1 showed that the macroscopic permeability of 
DMPC membranes for Co 2+ was altered by the insecticide lin- 



dane in agreement with the influence of the drug on the heat 
capacity profile. Other drugs as neurotransmitters also have the 
potential to shift melting transitions and to alter the fluctuations 
in membranes. This was for example shown for serotonine [57 1 
that has a quite significant influence on melting profiles. Thus, 
neurotransmitters also should display the potential to influence 
lipid channel events considerably. More so, every drug that al- 
ters the melting behavior of a membrane can change the appear- 
ance of lipid channels. This includes proteins and peptides that 
can alter lipid channel behavior without being channels them- 
selves, as discussed in section [8~2] 

The possibility to "block" or induce lipid channel events in 
protein-free membranes by drugs could be considered as an 
agonist or antagonists effect, or "channel-gating". It does not 
require specific binding to any macromolecule but rather repre- 
sents an unspecific physicochemical process. 

6.4 Voltage, pH and calcium 

About 10% of the lipids in most biomembranes carry negative 
charges in their head groups. Charged membranes display an 
electrostatic potential that depends on the charge density. At 
typical conditions of biological membranes, i.e., low charge 
density (10% of the lipids) and high ionic strength (wl50mM) 
the electrostatic surface potential of a membrane ^ can be ap- 
proximated by the low potential approximation that is given by 



1 a 

W = a oc —= 

e e K V c 



(32) 



where eo is the vacuum permittivity, e is the dielectric constant 
of water (w 80) and a = q/A is the charge density, k is the De- 
bye constant that is proportional to %fc with c being the ionic 
strength. Consequently, the larger the ionic strength and the 
smaller the surface charge density, the smaller the surface po- 
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tential. From this the electrostatic free energy can be calcu- 
lated: 



$> dq oc -—= 
„>=0 A Vc 



(33) 



i.e., the potential is quadratic in the total surface charge and in- 
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Figure 19: Currents through a membrane made of 
DOPC:DPPC=2:l mixture at 19° and 150mM KCl 
recorded at different voltages. One finds a linear current- 
voltage relationship for \V m \ <150 mV that becomes 
non-linear at higher voltages indicating that the voltage 
alters the state of the membrane. From H51V . 



versely proportional to the membrane area. Since biological 
membranes are asymmetrically charged (typically with most 
negative lipids being on the cytosolic side of the cell membrane 
[58 59 1) one expects a transmembrane voltage of order 50mV 
[60|. The molecular area of the gel state of a lipid membrane is 
about 25% smaller than that of a fluid lipid |48|. This implies 
that the electrostatic free energy of a gel membrane is larger 
than that of a fluid lipid membrane because the charge density 
is larger in the gel state. The decrease of the electrostatic free 
energy contribution upon membrane melting, AF e i l, leads to 



a lowering of the phase transition temperature by 

AF e , L 



AS 



(34) 



where AF e ix < is the electrostatic free energy difference 
between fluid and gel lipids. The mathematical expression 
for AF e i l (not given here) has been derived in [61 ] and been 
reviewed in lfl"3"l . The above equation implies that the change 
in lipid state from gel to fluid will lower the free energy of 
the membrane because the membrane area increases and the 
potential decreases. Further, by reversing the argument it is 
obvious that application of voltage across a membrane must 
change the lipid state. In particular, charged and uncharged 
membranes should respond differently to voltage. Since nega- 
tive lipid head groups can be protonated upon lowering the pH, 
the membrane state must be pH dependent. 
Voltage: The above argument showed that the lipid state 
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Figure 20: Effect of charge on the transition temperature 
of dimyristoyl methylphosphatidic acid membranes as a 
function of pH for four different ionic strengths (adapted 
from H611I ). The transition midpoint changes by about 14 
K towards higher temperature upon protonation of the 
membrane. The pKa of the membrane is at higher pH 
values when the salt concentration is lower. 



should depend on transmembrane voltage. An indication 
for this is shown in Fig{l9] that displays the transmembrane 
currents as a function of voltage, and the corresponding 
current-voltage relationship of the channels. It can be seen that 
the I-V curve is linear in a regime from about -150 to +150 mV 
(i.e., conductance is constant), while the I-V relation becomes 
non-linear at \V m \ > 150 mV. That voltage can reversibly in- 
duce channel events was reported by (29) ■ Antonov et al. |62| 
showed that voltage increases the melting transition in charged 
lipid membranes without significantly altering the transition 
temperature of zwitterionic lipids. However, in our daily lab 
routing we reproducibly find that zwitterionic membranes not 
displaying channels at low voltage show such events upon 
increasing the voltage (in agreement with eq. [TJ. Therefore, 
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we assume that voltage also has an effect on zwitterionic mem- 
branes. We further suspect that lipid ion channel formation is 
closely linked to the phenomenon of electroporation (cf. 11331 ). 
It should be added zwitterionic membranes are piezoelectric 
due to the electrical dipoles of their head groups [63 1. This was 
also shown in monolayer experiments demonstrating different 
dipole potentials in the gel and the fluid phases of the mem- 
branes and an effect of voltage on the lipid state |64|. 
The influence of voltage could be called 'voltage-sensitivity' 
or 'voltage-gating' . 
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Figure 21: Transmembrane currents through soy- 
bean phosphatidylcholine membranes in an unbuffered 
medium with 1M KCl. Traces are shown with different 
proton concentrations (pH range from 6.5 (trace a) to 1 
(traces f and g)) at different positive and negative volt- 
ages. Transmembrane currents are induced at low pH 
close to the pKa of the membrane. From MOV . 



Protons and pH: Let us consider a membrane that is 
negatively charged at neutral pH. Upon lowering the pH the 
phosphate groups will be protonated. This case has been 
treated in much detail by Trauble and collaborators |6T1 . They 
provided the electrostatic theory briefly outlined above and 
showed how the melting temperature depends on pH and ionic 
strength. Upon complete protonation the transition midpoints 
increases by 10-20 degrees. 



This is shown in Fig. 20 for 
dimyristoyl methylphosphatidic acid (DMMPA) membranes 
[61|. This is shown in Fig. [20] Similarly, zwitterionic lipids 
my become become positively charged at low pH values 
around 1-2. Their transition midpoint will decrease due to the 
electrostatic contribution of the lipid charges. From this it must 
be followed that the occurrence of lipid ion channels should 
be pH-sensitive. This has in fact been shown by Kaufmann 
and Silman ifTol . see Fig. 21 They found that channels can 



be induced by protons in the pH regime around 1 that is close 
to the pKa of the zwitterionic soybean phosphatidylcholine 
membrane. 



Calcium: Due to its two positive charges, calcium also dis- 
plays a pronounced effect on lipid membranes. It binds to both, 
negatively charged and somewhat weaker to zwitterionic mem- 
branes. This is likely to be associated with the electrostatics 
of the membranes. The two calcium charges probably cross- 
link the phosphate groups of two adjacent lipids. Typically, 
calcium increases melting temperatures. Depending on state of 
the membrane (e.g., whether found below or above a transition) 
calcium may therefore induce or inhibit lipid ion channels. This 
is shown in Fig. 
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Gogelein & Koepsell 11651 showed that 
channel events in brain phosphatidylserine can be "blocked" or 
inhibited by calcium (Fig. 22 top, while Antonov et al. |[T2l 



showed the reverse effect for DPPC membranes above the tran- 
sition (Fig. |22j bottom). The effect of calcium corresponds to 



OmV -50mV 



BmM CaClj 



2pA 



Is 




Figure 22: Influence of calcium on tlie occurrence of 
lipid channels. Top: Currents through bovine brain phos- 
phatidylserine : cholesterol=2:l (w/w) membranes in 
200mM cyclamate and pH 7.4. Addition of 8mM CaCl2 
inhibited the currents. From [65]. Bottom: Generation 
of quantized currents in DPPC membranes by calcium 
at T=64°C (the melting temperature of DPPC is 41°C). 
From HI 2V . 



a 'calcium-control' or 'calcium-gating' of the lipid channels. 
As for the case of drugs, no specific binding site at a particu- 
lar macromolecule is involved. One should bear in mind that 
both pH and calcium effects are related to the surface poten- 
tial and thus to the voltage effects. This includes the adaptation 
of lipid composition reported for many bacteria under various 
temperature, pressure or solvent conditions (see below). 

6.5 Other variables 

Everything that potentially changes the state of the lipid mem- 
brane will display the potential to generate or inhibit lipid ion 
channels. We just list a few of these: 
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• membrane curvature, because one of the work terms in 
the internal energy differential is related to curvature, and 
curvature has an influence on the membrane phase be- 
havior 1661. 



• dipolar interactions, because the work to orient a dipole 
(e.g., the lipid head group) in a field is part of the dif- 
ferential of the internal energy. It has been shown that 
electrostatic fields can alter the phase behavior of zwitte- 
rionic monolayers 11641 . 

• hydrolysis of membranes by enzymes, because hydroly- 
sis products change the phase behavior of the membranes 
||67l . Any change in membrane composition may have a 
similar effect. 



• membrane proteins. This case is discussed in section 8.2 



7 Channel lifetimes 

It has been found in many experiments that relaxation time 
scales are closely coupled to the melting process in membrane 
E3in]l22l|73l|31|75l|76l|57l.The fluctuation-dissipation the- 
orem reveals a fundamental coupling between the magnitude of 
fluctuations (expressed in heat capacity, compressibilities, ca- 
pacitance, ...) and the the relaxation time scales (indicating how 
fast a system relaxes to equilibrium after a perturbation), which 
are identical to the fluctuation time scales. Since we consider 
here the lipid channel formation process basically as the con- 
sequence of a local fluctuation in state, it is obvious that open 
and closed life times of a channel must be intimately related to 
the fluctuation time scale. Basically, the fluctuation-dissipation 
theorem states that at regimes of high compressibility and high 
heat capacity, fluctuations are slow. Thus, channel life times 
are expected to be long in the phase transition regime when 
their likelihood of formation is large. 



6.6 Free energy of channel formation 

Let us consider a membrane in a particular state that is not 
within the melting transition. This membrane displays a very 
small likelihood to display lipid ion channels. In order to in- 
crease the channel probability, the membrane has to be brought 
into the chain melting regime, which can be done be altering 
the various parameters described above. Heimburg & Jackson 
[ 53 1 showed that the free energy that is necessary to bring the 
membrane into a regime where channels occur (i.e., the free 
energy difference between the fluid and the gel phase) is given 
by: 



AG(T,x A ,Ap,...) 
( 

AH 



T m -T RT T 

rp XR XA+ryVAP T~ 

-L m. L\l± ± rn. 



\ v — 



(35) 



For each of the thermodynamic variables, there is one term in- 
fluencing the free energy difference. One could add terms for 
the electrostatic potential, for pH, calcium concentration, cur- 
vature, etc., which has not been done here. The content of this 
equation is: The further away a membrane is from its melting 
transition the more free energy has to be provided to generate 
channels. Further, for a given value of the free energy differ- 
ence the likelihood of channel formation is constant, but the 
values of the individual variables can change. The effect of 
anesthetics can be compensated by hydrostatic pressure, pH, 
or temperature change. These effects have been documented 
for anesthesia [68, 69, 53] yielding the well-known effects of 
pressure-reversal of anesthesia and the impossibility to anes- 
thetize inflamed tissue with normal anesthetic doses. 



7.1 Theoretical considerations 

We now return to the theoretical treatment of fluctuations given 
in sections [572| and [531 On the basis of this description the fluc- 
tuation life times of membranes can be obtained |76ll57l . Lars 
Onsager [44 1 treated time dependent phenomena such that he 
assumed a proportional relation between fluxes (rates of change 
of a variable) with the thermodynamic forces, an assumption 
corresponding to the equation of motion in viscous fluids where 
the influence of inertia is small. Following this approach, the 
flux of heat (enthalpy) of the membrane after a perturbation 
back to equilibrium (which could consist of the thermal colli- 
sions related to Brownian motion) is given by 



dt 



L- X 



H 



(36) 



where Xh is the thermodynamic force related to enthalpy dif- 
ferences (cf. section 5.2 1, given by DS/dH with the thermo- 

H 



section 
dynamic variable 



= H — (H). L is a material constant 
that we assume to be temperature independent. This equation 
is a simple version of Onsager's phenomenological equations if 
only one independent thermodynamic variable is present. 



X 



H 



dS 
dH 



-g HH -(H-(H)) 



According to eqs. 10 12 and 13 
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This leads to 



d(H-(H)) 
dt 



HH 



L 

TV 



or g HH = 



(H-(H)) 



T 2 c 



(37) 



(38) 



(39) 
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with the solution 

(H-(H))(t) 



(H - (H)) ■ exp 



T 2 c 



(H - (H)) ■ exp - 



introducing the relaxation time 



L 



(40) 



(41) 



This equation implies that the relaxation (or fluctuation) time 
scale is large when the heat capacity is high. If we consider 
the channel formation as a fluctuation, we have to conclude 
that the typical time scale of channel opening and closing is 
proportional to the heat capacity. In particular, in the transition 
where the likelihood of lipid ion channel formation is high, they 
should also display long open times. 

It has been found that for pure lipid membranes the phe- 
nomenological coefficient L is of order 6-20- 10 8 J K/mol sec. 
For a DMPC multilayer one finds a relaxation time of about 30 
seconds at the maximum. Transitions in biological membranes 
are much broader and the heat capacity maximum displays 
much smaller values. Assuming the same phenomenological 
coefficient as for DMPC membranes, one estimates about 100 
milliseconds at the c p maximum of bovine lung surfactant or 
E.coli membranes. Since biological membranes display a c p 
maximum some degrees below body temperature, the expected 
open times at physiological temperature is expected to be of 
the order of a few milliseconds. Interestingly, this time scale is 
just of the same order as the dwell times of quantized current 
events in section |4] 

Summarizing: While the channel open times in artificial 
membranes at the c p maximum may be very long (of order sec- 
onds), they are expected to be of millisecond order for biomem- 
branes at physiological temperature. This time scale, expected 
for lipid channels, is interestingly also exactly that found for 
protein channels. This is unlikely to be accidental. 

7.2 Experiments 

It has been shown by 1 30 ] that channel life times are in fact 
much longer in the transition of a lipid mixture that displays 
a transition half width comparable to that of biomembranes. 



This is shown in Fig 23 Both heat capacity profile and mean 
conductance of a lipid mixture of DC15PC and DOPC (95:5) 



display maxima at about 30°C (Fig 23 top). The typical open 
time of the lipid channels is of order 20 milliseconds, while it 
is rather 3ms above the transition. If one compares the time 
scales of the current fluctuations in the phase transition regime 
and in the fluid phase one finds that they are about one order 



of magnitude different (Fig 23 bottom), in approximate agree- 
ment with the above calculation (eq. 41 1. The timescale in 



1.2 
1.0 
08 
0.6 
04 
0.2 



phase transition 




fluid phase 




. a gel phase 






20 30 40 




Temperature [ a C] 








M 



20 25 30 35 
Temperature [°C] 



40 



45 




Figure 23: Top: Relative conductance of a D15PC- 
DOPC (95/5) mixture as a function of temperature at V m 
= 90mV. The insert shows the melting profile of the lipid 
mixture. Adapted from H30V . Bottom: Logarithmic plot of 
the timescales of the current fluctuation in the fluid phase 
(circles) and the phase transition regime (triangles). The 
average lifetimes center around 3ms in the fluid and 20ms 
in the phase transition, which is in good agreement with 
our theoretical prediction. From H301I . 



transition it can only be approximated by a double-exponential 
function. Gallaher et al. 117711 provided an analysis showing 
that the current fluctuation time scales in the transition regime 
obey a power law rather than single exponential behavior as 
suggested by the derivation above. The authors rationalized 
this by assuming that the conduction pores may exist in dif- 
ferent states with different gel-fluid environments. Figs. 13 



the fluid phase follows a single-exponential decay while in the 



and 14 also suggest that there might be more than one possible 
pore geometry with similar pore dimension. 
Summarizing, it seems as if one can understand the time scale 
of lipid ion channel opening and closing with simple non- 
equilibrium thermodynamics (fluctuation-dissipation). This 
means that heat capacity, phase behavior and time scales are 
different aspects of the same physics. Whenever heat capacities 
and permeabilities are high, time scales are larger. However, 
while the overall change in time scales is consistent with the 
above derivation, multi-exponential or power law behavior is 
not contained and requires more sophisticated analysis. 
Heat capacity profiles can be influenced by anesthetics, neu- 
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rotransmitters and by membrane peptides. Interestingly, it has 
been shown by 071 that those c p changes are accompanied 
by equivalent changes in the relaxation time scales, such that 
those molecules cause a predictable change of the relaxation 
time scales. The role of proteins is outlined in more detail in 
the following section. 



8 Interfaces: Domains and proteins 

Fluctuations are large in regions where the free energy of gel 
and fluid lipids are similar. Therefore, they are maximum at the 
melting transition. As shown above, many variables can be ad- 
justed (pressure, temperature, ...) such that the membrane sys- 
tem is situated in this regime. However - such fluctuations must 
not be the same everywhere in the membrane if the membrane 
is heterogenous. Biological membranes are complex mixtures 
of lipids with different melting points and different hydropho- 
bic lengths, and many other molecules including cholesterol 
and proteins. These membranes display phase separation and 
domain formation phenomena. Therefore, one must expect that 
also the fluctuations display different magnitude in different 
membrane regions. It is especially interesting to consider in- 
terfaces between fluid and gel domains. 

The magnitude of fluctuations depends on the length scale over 
which they are monitored. One can consider fluctuations of 
the whole membrane, which enter into the heat capacity or the 
macroscopic compressibility. If one considers the fluctuations 
in windows of smaller size, e.g. on the length scale of indi- 
vidual domains or on molecular scale, one obtains local heat 
capacities and compressibilities. There are different ways of 
probing the local fluctuations, e.g., 

• experimentally: One can investigate the softness of a 
membrane by AFM. One obtains the compressibility on 
the length scale of an AFM-tip (around lOnm). This has 
been done for peptide containing membranes in |78 1. 

• theoretically by simulations: In Monte Carlo simulation, 
one can directly obtain the fluctuation on arbitrary scales. 
Such simulations typically use experimental observables 
such as melting enthalpy and entropy, and the shape of 
the melting profile in order to determine the simulation 
parameters (e.g., |[79] [80] [8T] [82l |83]|). In principle, such 
information is also contained in molecular dynamics 
(MD) simulations. However, to our knowledge MD has 
not been used yet to obtain such information. 



8.1 Domain interfaces 



In Fig. 24 domain formation of lipid mixtures is shown. Confo- 
cal microscopy images of large unilamellar vesicles display do- 
mains of gel and fluid state (Fig. [24] top). Such domain forma- 
tion can be described by lattice Monte Carlo simulations (Fig. 
24] bottom left). Details are, e.g., given in l85l . For the purpose 
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Figure 24: Top: Confocal microscopy image of a DLPC- 
DPPC (1:2) mixture in its melting regime at 27° C. Dark 
regions represent gel domains, bright regions are fluid 
H84V . Bottom, left: Snapshot of a Monte Carlo simulation 
of a DMPC-DSPC (1:1) mixture in the melting regime 
(36.9° C). Dark shades correspond to gel domains while 
bright shades represent fluid domains. Bottom, right: Lo- 
cal fluctuations of the same matrix. Bright shades corre- 
spond to large fluctuations while dark shades correspond 
to small fluctuations. The magnitude of the fluctuations 
is related to the membrane permeability and the likeli- 
hood of lipid channel formation. For the example shown, 
fluctuations are larger in the fluid regions than in the gel 
regions, but are maximum at the interface between the do- 
mains. These regions are those of maximum permeability. 
From l[85\l . 



of this review only the results are important. The fluctuations 
of single sites in this matrix (local fluctuations) are shown in 
Fig. 



24 (bottom right). It is found that under the conditions 



indicated in the figure legend the gel domains display much 
smaller local fluctuations than the fluid domains. The largest 
fluctuations, however, are seen at the domain interface. Thus, 
while the domain interfaces are the regions of largest perme- 
ability, also the domains themselves contribute to permeation 
(in Fig|24] (bottom right) especially the fluid domains). This 
sheds further light on the two possible explanations for the lo- 



cation of permeation pores discussed in section 5.5 
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8.2 Proteins as membrane perturbations 

The gel-fluid domain interface is not the only region of large 
fluctuations. Macromolecules as proteins can also create in- 
terfaces with lipids that influence their phase behavior. Two 
examples are given in Fig. 25 showing the influence of cy- 



tochrome bs and band 3 protein of erythrocytes on the melting 
behavior of DPPC and DMPC membranes. While cyt. b$ 
lowers the melting point, band 3 protein is increasing it. It 
has been shown in [83| that this can be understood if one 
assumes that the first protein mixes well with fluid membranes 
but not with gel membranes, while the second protein mixes 
better with gel and not so well with fluid regions. Mouritsen & 
Bloom [ 88 1 proposed that this is due to hydrophobic matching 
that compares the length of the hydrophobic core of the integral 
protein with the thickness of the membrane. Good miscibility 
suggests that the hydrophobic core of the protein has a similar 
length then the surrounding membrane. Bad miscibility indi- 
cates hydrophobic mismatch. 

Gramicidin A is a small dimeric channel peptide. As shown 
in Fig. 26 (top) it is significantly shorter than a gel lipid of 



DPPC. Thus, one expects unfavorable hydrophobic matching 



in the gel phase of DPPC. Fig. 26 (botto left) shows the distri- 
bution of gel and fluid state lipids and the proteins as obtained 
from a Monte Carlo simulation that was based on heat capac- 
ity analysis J83). The simulation was performed slightly below 
the melting temperature of DPPC. While most of the matrix is 
found in its gel state with some remaining fluid domains, the 



proteins aggregate into smaller clusters. Fig. 26 (bottom right) 
shows the corresponding local fluctuations. One finds that the 
fluctuations are maximum at the interface of the protein clus- 
ters. This indicates that the protein interface also forms a mem- 
brane region with larger likelihood of lipid channel formation. 



The results from the simulation are summarized in Fig. 27 that 
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Figure 25: Influence of proteins on the melting behav- 
ior of lipid membranes. Left: Cytochrome b$ lowers the 
melting point (from i T73l/ adapted from j|56l/). Right: Band 
3 protein of erythrocytes increases the melting point (from 
Kl3]l adapted from [87]). 



is a variation of a picture from Mouritsen & Bloom used to 
describe their mattress model [88 1. If a protein with large hy- 
drophobic length is penetrating a fluid membrane with smaller 
thickness, the lipids at the interface of the protein will be close 
to the melting regime, which is due to wetting of the hydropho- 



bic core of the protein (Fig. 27 top). Thus, at this interface fluc- 



tuations are large and the likelihood of lipid channel formation 
is maximum. Similarly, if the protein has a short hydrophobic 
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Figure 26: Simulation of tlie antibiotic peptide grami- 
cidin A in a DPPC membrane. Top: Gramicidin A is 
very short as compared to the gel lipids of DPPC and 
displays an unfavorable hydrophobic matching. Bottom, 
left: Monte Carlo snapshot of a simulation of gramicidin 
A in a DPPC membrane. Slightly below the melting tran- 
sition one finds the membrane mostly in the gel state (gel 
lipids - dark grey) with some remaining fluid domains 
present (fluid lipids - light grey). Due to the unfavorable 
hydrophobic matching, the proteins aggregate in the gel 
phase into clusters (black dots - gramicidin molecules). 
Bottom, right: For the same snapshot the magnitude of 
the fluctuations is shown locally. Dark grey shades corre- 
spond to small fluctuations while bright shades represent 
large fluctuations with a high likelihood of pore forma- 
tion. 
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Figure 27: Schematic drawing of proteins as perturbation 
of the lipid membrane state. Top: Protein with long hy- 
drophobic length embedded into a fluid lipid membrane. 
Bottom: Protein with short hydrophobic length embedded 
into a gel lipid membrane. Due to wetting of the protein 
interface the lipids close to the protein are in a transition 
state with high fluctuations and an increased likelihood 
of lipid channel formation. 



length and is surrounded by a gel membrane of larger thickness, 
the protein interface will constitute a location for lipid channel 
formation (Fig. 27 bottom). It is a general predictions from 



such simulations that the tendency of proteins to aggregate and 
their potential to trigger lipid ion channel formation are closely 
coupled. 

The role of the protein in such a scenario can be seen as dop- 
ing the membrane and influence the physics of its surrounding, 
much like in semiconductors where the doping agent gallium 
can influence the physics of silicon. However, the proteins in 
this situation do not conduct ions themselves. Rather is the 
outer interface of the protein the site of permeation. All argu- 
ments arguing in favor of a selectivit^of a protein pore (chan- 
nel protein) would also be true for the outer interface. Thus, 
the protein in the above situation could be considered as an in- 
verted pore (the pore interface is on the outside). The physics, 
however, is that of the lipid membrane. 

9 Discussion 

We have reviewed here some of the literature indicating the 
presence of quantized channel events in lipid membranes. It is 



a striking fact that one finds such events in membranes that do 
not contain proteins. The currents display typical amplitudes in 
the lOpA regime at voltages around lOOmV. Typical life times 
are of order 10 milliseconds but can be much longer in some 
experiments. These currents can be inhibited by anesthetics [2 1 
but most likely by everything that shifts melting transitions, in- 
cluding peptides and neurotransmitters. They respond to pres- 
sure, tension, pH, calcium, voltage and other variables. In some 
sense one could consider the overall membrane as a sensitive 
receptor for all of these variables if one is close to a transition. 

9.1 Comparison with proteins 

When comparing the lipid ion channels with protein channels 
it has to be stated that 

• typical conductances and lifetimes of protein and lipid 
channels are of similar magnitude. The two classes of 
events are seemingly indistinguishable. 

If two events cannot be distinguished one has to reconsider 
what the factual evidence actually is that allows to attribute 
the current events to different processes. The findings on lipid 
channels reported here imply that 

1 . the observation of quantized currents is no proof for the 
action of ion channel proteins. 

2. the action of a drug on a conduction event is no proof for 
specific binding of the drug to a channel protein. 

In logics, elementary rules of inference distinguish between 
necessary and sufficient criteria. For channel proteins it must 
be concluded that the finding of quantized currents alone is 
neither necessary nor sufficient to infer the action of a receptor. 
The quantized current could be though the lipid membrane, 
and there is no particular reason why the a quantization should 
indicate the action of a macromolecule. Similarly, the action of 
drugs is not sufficient to infer specific binding. The drug could 



3 According to Einstein 1 89 1 an all-or-nothing selectivity of a membrane (or 
a membrane protein) for a particular substance violates the laws of thermody- 
namics. Therefore, we rather refer to different permeation rates. 



have an effect without binding, as shown in section 6.3 for the 
anesthetics. For example, the inference that the abolition of 
channel events in biomembranes by tetrodotoxin (TTX) proves 
the action of a sodium channel protein and the specific binding 
of TTX to this channel is not admissible without additional 
information, which is typically not provided in the literature. 

The experiments reviewed in this paper do not give an in- 
dication for that protein channels do not do exactly what has 
been suggested in the protein ion channel literature. However, 
the data provide considerable evidence that one can find the 
same events in the absence of proteins and therefore one does 
not necessarily require proteins as an explanation. While the 
existence of lipid channels is easy to prove because a lipid 
membrane can be measured in the absence of proteins, it is 
presently impossible to measure protein channel conductances 
in the absence of a membrane. Thus, it is much more difficult 
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to actually prove that a protein ion channel is really active in an 
experiment. From a patch-clamp experiment it is impossible to 
conclude, which path a particular current event has taken. The 
patch area is between 1-100 fim 2 while the cross section of a 
single channel protein, e.g. the potassium channel, is rather 25 
nm 2 = 25 • 10 _6 ^m 2 . A typical aqueous pore has a diameter of 
0.5 nm. Since no experiment can actually monitor the path of 
the ions and therefore demonstrate the existence of a particular 
pore inside a membrane protein, the argument is typically of 
rather indirect nature: 

• the protein increases the probability to find conduction 
events, and conductivities are small in the absence of the 
proteins 

• the currents can be blocked by specific drugs as 
tetrodotoxin (TTX) or tetraethylamonium (TEA) 

• site specific mutation can alter or block the currents 

One must state that it has been incorrectly assumed that the 
lipid membrane is inert and does not display quantized con- 
duction events. Second, proteins can act as perturbations of the 
membrane and thus have the potential of altering the state of 
the membrane. A direct proof for this effect is still missing, but 
it is straight-forward to conclude this from the thermodynamics 
described in this article. Therefore, a conductance that displays 
a protein concentration dependence must be taken as a proof 
for a correlation of the proteins with the conduction events but 
not as a proof for the conductance through a pore in the protein 
itself. Conduction events in lipid membranes depend also on 
the concentration of drugs like anesthetics (as shown in section 
6.3 1 even though anesthetic molecules do not form channels. 
The dangers of confusing lipid with protein channels may be 
evident in a publication by Cannon et al. |90|. These authors 
found that calcium channels reconstituted in lipid mixtures are 
not only most active at the phase boundaries of the lipid phase 
diagram but also display the longest dwell times. While there 
is no particular reason why a protein should behave this way, 
this is exactly the behavior expected for lipid ion channels in 
lipid membranes at phase boundaries (dwell times are long 
in the transition, see section [7J. The channel conductances in 
this paper were the same as, e.g., in the lipid mixture shown in 
Fig. 
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Further, MacKinnon and collaborators |91| showed 
that the potassium channel conduction properties depend very 
sensitively on the lipid environment proving a connection 
between state of the lipid membrane and protein channels. The 
authors assumed that the voltage sensors of these proteins are 
lipid-controlled. The effect of anesthetics on lipid membrane 
pores is very similar to that reported for the acetylcholine 
receptor and many other channel proteins (reviewed in ifBTI . 



9.2 Nerve pulses 

In a series of recent publications it has been shown that close 
to melting transitions of lipid membranes one expects the 



possibility of electromechanical pulse propagation (solitons) in 
cylindrical membrane tubes [20, 92 60 93 1. These solitons in 
fact share many properties with the action potential in nerves 
that cannot be explained by the Hodgkin-Huxley model, e.g., 
reversible heating of the nerve, and mechanical alterations like 
thickness changes and length changes. Thus, it may be pos- 
sible that the nerve pulse can be explained without involving 
ion channel proteins. However, exactly under the conditions 
under which such electromechanical solitons can propagate 
one also expects the formation of lipid ion channels. As shown 
in this article, everything that alters the the phase behavior will 
alter the lipid channel appearance and simultaneously the free 
energy necessary to induce a nervous pulse when following an 
electromechanical approach to rationalize the phenomenology 
of the action potential [20|. This is a beautiful correlation 
purely based on the thermodynamics of the lipid membrane 
that is well understood as shown here. 

9.3 Conclusion 

Taking the extreme relevance of the finding of quantized cur- 
rents in artificial membrane where any effect of proteins can 
be excluded it is quite surprising that this phenomenon has 
not been studied more intensively. Not only do these currents 
display similar conductances as proteins channels, but obvi- 
ously similar fluctuation life times. Since many biomembranes 
exist in a state close to melting transitions of their membranes, 
the striking similarity between the lipid and protein events 
suggests that these events are difficult or even impossible 
to distinguish. Should in fact some of the reported protein 
conductances be due to lipid pores, the theoretical description 
of such channels would rather lie in the thermodynamics of 
the membrane and its cooperative phase behavior than in the 
geometry of individual proteins. A discussion of the lipid 
channel phenomenon and its connection with protein channels 
is overdue. 

Even though there have been reports on lipid ion channels 
for 30 years the investigation of lipid ion channels is still 
in an early stage. Many of the thermodynamics couplings 
described in this article have not been tested yet by experiment. 
However, already now many of the implications of applying 
the fluctuation-dissipation theorem to lipid channels found 
experimental verification, e.g., the coupling between life times 
and heat capacity, the effect of pH, anesthetic drugs and lateral 
tension. Many further experiments are needed to get better 
statistics of the membrane conductance and the influence of 
various parameters. It would in particular be nice to demon- 
strate the influence of proteins that cannot be channel proteins 
(for instance peripheral proteins) on the quantized conductance 
of lipid membranes. Such experiments are not difficult and 
will most likely be done in the near future. 
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